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Abstract 



We consider a general symmetric (3 x 3) mass matrix for three generations 
of neutrinos. Imposing the constraints, from the atmospheric neutrino and 
solar neutrino anomalies as well as from the CHOOZ experiment, on the mass 
squared differences and on the mixing angles, we identify the ranges of allowed 
inputs for the 6 matrix elements. We apply our results to Majorana left- 
handed neutrino masses generated at tree level and through fermion-sfermion 
loop diagrams in the MSSM with R-parity violation. The present experimen- 
tal results on neutrinos from laboratories, cosmology and astrophysics are 
implemented to either put bounds on trilinear (Aj^fc, A'^^-^) and bilinear (/ie./^.r) 
R-parity-violating couplings or constrain combinations of products of these 
couplings. 
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I. INTRODUCTION 



The current experimental evidence strongly suggests the existence of a non-trivial struc- 
ture of the mass matrix for 3 generations of neutrinos. Indeed, deficits of the solar electron 
neutrino and atmospheric muon neutrino fluxes compared to their theoretical predictions 
are strong evidence in favour of neutrino oscillations as an explanation of these deficits. 
Furthermore, oscillations of neutrinos imply masses for neutrinos and mixing angles, which 
relate the flavour basis to the mass eigenstate basis, as is the case in the quark sector. 

Observations in favour of the oscillation solutions came recently with the Super- 
Kamiokande collaboration results, which have to be added to those of other atmospheric 
neutrino experiments (1MB [Q], Soudan 0, Kamiokande [^) and solar neutrino experiments 
(Homestake §, Gallex §, SAGE [0], Kamiokande [§, Super-Kamiokande §, MACRO 
and LSND [|ll|). The oscillation explanation [| of the solar neutrino problem, the atmo- 
spheric neutrino anomaly and the LSND results suggest three very different values of neu- 
trino mass squared differences, namely Am^^j^ ^ Am^^^ <^ AmLgNO) "with Am^^^ < 10~^ 
eV^, Aml,^ E [10-3, IQ-^] eV^ and Am^gNo e [0.3, 10] eV^. The evidence for 
observed by LSND has not been confirmed or excluded by the KARMEN [J6[ experiment. 
The MiniBooNE at FNAL or MINOS long-baseline experiments could provide the answer. 
Given the atmospheric and solar neutrino observations, and excluding the LSND results 0, 
the deduced pattern of neutrino masses and mixing requires two different neutrino mass 
squared differences and three mixing angles. Recent interest in this subject has produced 
a very large number of possible models, which can accommodate the data, see refs. for 
more details. 

Our approach differs from most of the analyses in that we do not specify the model 
initially. In a general manner, we will use the current relevant information on the mixing 
angles and mass squared differences to constrain the neutrino mass matrix elements in the 
flavour basis (z/g, z/^, i^r). Consequently this will constrain the parameters of the involved 
model. The results of this work can be applied to any class of model with three neutrino 



^ Other possible solutions, which could accommodate the experimental neutrino data, can involve 
flavour-violating interactions, or other more exotic possibilities such as the violation of Lorentz 
invariance [^ or the equivalence principle [13|. However, in the case of atmospheric neutrinos, it 
has been shown [14|, [15| that neutrino decays or modification of special or general relativity do 
not fit the data. 

^In a forthcoming analysis, we will include LSND results, which calls for a fourth neutrino species 
or a sterile neutrino. 
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generations. We also present an application of our results of the specific Minimal Supersym- 
metric Standard Model (MSSM ) with R-parity violation. We will also try to elucidate if any 
of the different possible neutrino spectra (degenerate neutrinos, hierarchical, pseudo-Dirac 
p!8| ) are preferred for each solution of the solar anomaly combined with SuperKamiokande 
and CHOOZ constraints. Even barring the presence of a sterile neutrino, other interesting 
bounds on neutrino masses arise from cosmological or astrophysical sources, such as hav- 
ing a neutrino component of hot dark matter (HDM). In addition, low energy experiments 
can also constrain neutrino masses; the strongest constraint, which is applicable when the 
neutrinos are Majorana particles, arises from neutrinoless double beta decay. We will also 
analyse our results in view of the two latter bounds. 

The neutrino is massless in the Standard Model and, if massive, requires the presence of 
new physics. For reviews, concerning models for neutrino masses and neutrino properties, 
see [p!9|-p^ and references therein. If the neutrino has a Majorana mass, the model must 
violate lepton number conservation. In supersymmetric (SUSY) extensions, gauge invariance 
and renormalizability no longer ensure lepton number L (or baryon number B) conservation. 
The generalization of the MSSM, which allows R-parity (Rp) violation (Rp = (_i)^+3B+2S^ 
where L, S are the lepton and baryon number and the spin of the particle, respectively), 
allows left-handed neutrinos to obtain a Majorana mass, at tree level through mixing with 
the neutralinos, and through loop diagrams that violate lepton number (in two units). We 
apply our general results to the MSSM with R-parity violation at one-loop order, allowing the 
presence of both bilinear and trilinear Rp- violating couplings without imposing any hierarchy 
between the parameters. We thus constrain all relevant Rp-conserving and non-conserving 
couplings by imposing the constraints from atmospheric and solar neutrino anomalies and 
the CHOOZ experiment. As a by-product, we see that in this model the constraints from 
neutrinoless double beta decay and cosmology are automatically satisfied. 

The paper is organized as follows. In section II we expand on the inputs used. In section 
III we establish our notation for the general case of a 3-generation neutrino mass matrix. 
In section IV we introduce our notation for an application of our general results to the 
MSSM with Rp-violation. We also give the results for the tree-level neutrino mass matrix 
and the loop corrections to each matrix element. We start by parametrizing the general 
neutrino eigenmasses and eigenstates and use the neutrino data to constrain the ranges of 
variation of the general inputs. In section IV.A we introduce a toy-model that corresponds to 
several limiting cases of the MSSM without R-parity conservation. In section V we present 
our results for the general case and for the specific case of Majorana masses arising in the 
R-parity violating model. Finally in section VI we summarize and conclude. 
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II. GENERAL INPUTS 



In this work, we compute the masses and mixing angles of the neutrino spectrum assum- 
ing three generations of neutrinos (this means that we are not taking the LSND result into 
account). We then apply all relevant existing bounds on the neutrino masses and mixing 
angles to constrain the elements of our neutrino mass matrix in the flavour basis. These 



inputs are from laboratories such as the CHOOZ experiment [^], the relevant mass squared 
differences Am^ and mixing angles sin^ 26 assuming the oscillation solution as an explana- 
tion of the solar and atmospheric neutrino deficits. We then check whether the upper bound 
coming from neutrinoless beta decay {f3P)ou [p^-p6| and also the global cosmological upper 
bound on the sum of the neutrino masses |^ are satisfied. The best limit on the effective 
mass nies appearing in (/3/3)oy, and defined by 



rrics 



(1) 



has been derived in the Heidelberg-Moscow ^^Ge experiment [^: |meff| < (0.5 — 1.0) eV 



(at 90% C.L.). The CP phases that might appear in the mixing matrix Uek are not relevant 
in our derivation, because we only use the second r.h.s. term in eq. (|l|). Another direct upper 
bound comes from the measurement of the high energy part of the P-decaj spectrum: the 
upper limits on the electron neutrino mass obtained in the TROITSK and MAINZ p9[ 
experiments are m,^ < 2.5 eV and m,^ < 3.4 eV, respectively. However, these experiments 
suffer from some ambiguities referred to as "the negative mass squared problem" , which is 
still not completely understood and we will therefore not use their bounds. 

There are also upper bounds on the neutrino masses coming from astrophysical and 
cosmological considerations such as the one from hot dark matter (HDM), which suggests 
that J2i = mi + 1712 + m-s < few eV 0. 

Global fits to the data show that neutrino oscillations among three neutrino flavours 

p, n, m, 11 



are sufficient to accommodate the solar and atmospheric data |1 
Three different types of solar experiments (using different detectors: chlorine, gallium and 
Cherenkov) are sensitive to different solar neutrino energy ranges; to accommodate them, 
three regions of oscillation parameter space are allowed [31|, |34|, [BH], which correspond to 



the vacuum oscillation solution, MSW with a large mixing angle (MSW-LMA) and MSW 



recent analysis (observation of distant objects favouring a large cosmological constant instead 



of HDM) suggesting a non-zero cosmological constant |3C] shows that the bound of a few eV for 
the sum of the neutrino masses is no longer imposed, since a hot dark matter component is not a 
necessary ingredient in this case. 
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with a small mixing angle (MSW-SMA). The estimates required to fit the data for the neu- 
trino mass squared differences and the mixing angles are shown in table I. It is important to 
keep in mind that the neutrino oscillation scenario can only restrict mass squared differences 
but not fix the absolute neutrino mass scale; moreover, it cannot distinguish whether the 
smallest mass splitting, which corresponds to the solar data, is between the two lightest 
mass eigenstates or the two heavier ones. 

To fix our notation, for this case with three generations, the flavour states ui are expressed 
in terms of the mass eigenstates z/j using the 3x3 mixing matrix U 

y, = j^Ul^^. (2) 

i=\ 

We hierarchically order the mass eigenvalues and denote them rrij, z = 1,2,3, such that 

mi < 777,2 < "'^3- 

In view of the ranges of the mass squared in table I, if the neutrino spectrum is hierarchical, 
then A?77^2 = vr?^ — is the parameter relevant to the solar case and Am^g = mg — vn\ 
is the dominant one for the atmospheric case with Aml2 <^ Am^g and Am^g ~ Am^g. 
The corresponding mixing angles (see the parametrization below) are sin 2^12 and sin^ 2^23, 
respectively. Besides the hierarchical mass spectrum, there exist two other alternative neu- 
trino spectra [T^, one in which the three neutrinos are quasi-degenerate and another 
(called pseudo-Dirac) in which only two neutrinos are quasi-degenerate and the third is much 
lighter (or heavier). In table II we show the different possible regimes and indicate which is 
the corresponding mass squared difference. 



Experiment 


Am^ (eV^) 


sin^ 26 


Atmospheric 


(1 - 8) X 10-3 


0.85 - 1 


Solar 






-MSW-LMA 


(3 - 30) X 10-5 


0.6 - 1 


-MSW-SMA 


(0.4-1) X 10-5 


10-3 - 10-2 


-Vacuum 


(0.5 - 8) X 10-1° 


0.5 - 1 


CHOOZ 


> 3 X 10-3 


sin 6*13 < 0.22 



TABLE I. MSW-LMA, MSW-SMA and Vacuum stand for MSW large mixing angle, small 
mixing and vacuum oscillation solutions, respectively. 



III. COMPUTATION OF THE MASS AND MIXING MATRIX 
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Spectrum 


Solar 


Atmospheric 


Hierarchy 




Amfg 


Degenerate 


Am23 or Amfg 


Amf3 


Pseudo-Dirac 


Am23 ^^^12 


Amf3 



TABLE II. Different possible regimes and corresponding mass squared difference. 



A. Generic mass matrix 



We consider a general 3x3 symmetric mass matrix in the flavour basis for three gener- 
ations of neutrinos (z/g, ly^, Vr)'- 



M 



Mi2 M22 M23 
[Mi3 M23 M33y 



(3) 



We construct the mass squared differences of the eigenvalues of Ai and the mixing matrix, 
which relates the flavour basis to the mass eigenstates, in a general way in terms of these 
input elements. We denote the three mass eigenvalues by rrii for i = 1, 2, 3, which satisfies, as 
mentioned above, the hierarchy mi < m2 < m^. For any given model, these masses and also 
the mixing matrix would be given in terms of the parameters of the model. We then apply 
the constraints, given in table I, individually and combined. The combined constraints are 
obtained by requiring the atmospheric, CHOOZ and one of the solar solutions to be satisfied 
simultaneously. In this way we can determine the region, in the (Am^ , sin^ 29) plane, that 
satisifies each of the three possible combinations of constraints simultaneously. The answer 
is then translated into terms of ranges where the parameters of the model vary so as to fulfil 
neutrino data constraints. 



B. Generic mixing matrix 



Computing the mass eigenstates of the mass matrix (|^), we obtain the rotation matrix, 
which contains the relevant mixing angles. We use the Chau and Keung parametrization of 



a 3 X 3 matrix |37| given by 



U 



C12 Cl3 

-S12 C23 - C12 S23 Si3e 



Sl2 Ci3 

C12 C23 - S12 S23 sise 



1(5 



Sue-'' \ 

S23 Ci3 



diag{e^"^e^°M}, (4) 



V Si2 S23 — C12 C23 Si3e* — C12 S23 — S12 C23 Si3e* C23 Ci3 / 



where c^- 



cos(6'jj) and Sij 



sm 



5 is the Dirac CP phase and ai^2 are the Majorana 
ones (we have two additional CP phases in the case of Majorana particles). As previously 
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mentioned, we are going to work with a real mixing matrix and will ignore CP phases because 
we are interested in our derivation in the modulus of the matrix elements, as we will see. 

The survival probability Pu^^u^y relevant to the case of solar fluxes, depends only on the 
first row of the mixing matrix in eq. (^), i.e. on |f/ejp, with i = 1,2,3. In the atmospheric 
case, the oscillation probability depends on the last column of (^, i.e. on |t/^3p, with 
£ = e, /i, T. The other elements of the matrix are not constrained by any direct experimental 
observation. With this parametrization, we directly obtain the parameters sin^ 29i2 and 
sin^ 2^23 and also the relevant CHOOZ parameter sin ^13, which we need for our general 
analysis. 



IV. APPLICATION: MSSM WITH R-PARITY VIOLATION 



The specific Rp-violating model to which we apply our general results has been exten- 
sively analysed in the context of pure phenomenology, as well as in trying to fit the data of 

H, H, 0, glj, H, [H, 0, H, 0, 13 different 



neutrino anomalies. In refs. 
phenomenological aspects of this theory have been analysed with respect to neutrino masses. 
In particular, refs. pO| , ^ have established in a basis-invariant way that neutrino masses are 
always generated in these models even when universality of the soft SUSY-breaking terms 
is assumed at some high scale. Previous studies have also tried to constrain the different 
Rp-violating couplings that appear in the MSSM Lagrangian, considering only the effect of 
bilinear terms |51[| or only of trilinear couplings [|5^, or from both |6^, solar 
and atmospheric neutrino data. Both tree-level and one-loop effects have been considered. 

As stated in the introduction here we do not bias our selection of Rp-violating couplings, 
which would be contrary to a common practice of keeping only one R-parity-violating pa- 
rameter at a time. We also include all possible bilinear and trilinear parameters for the 
MSSM without R-parity. In addition we do not impose any hierarchy between the trilinear 
couplings. Finally, our results present the allowed ranges for these parameters using the 
solar and atmospheric anomalies and the result from the CHOOZ experiment only. We will 
present our comparison to previous analyses in the last section. 

The most general renormalizable superpotential for the supersymmetric Standard Model 
with lepton-number violation is 



k J) 



(5) 



the {i,j,k) are flavour indices, (a, 6) are SU{2) indices, and the {a, (3) are flavour indices 
running from to 3. The are the doublet superfields with hypercharge Y = —1. Note that 
the A couplings are antisymmetric in the first two indices. The usual R-parity-preserving 
Lagrangian is obtained when only /Xo, Xoik = Kok = ^tk non-zero and we can identify 
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Lo = Hd- In the model of eq. (|^) we have 9 additional A couplings and 27 new A' couplings 
compared to the R-parity-conserving case. Note that thanks to the additional degrees of 
freedom, we can rotate in the flavour space of the "down-type" scalar fields to set the vacuum 
expectation values of the sleptons to be zero 0. Henceforth, we will work consistently in this 
basis. 

We use the following parametrization for the tree-level mass matrix in the basis of 
charged-lepton eigenstates: 



tree 



2 

d^ dfi^T 



(6) 



the corresponding expressions for the ctj's are given in appendix]^. This tree- level mass ma- 
trix has only one non-zero eigenvalue. It is thus necessary to include additional corrections, 
which could provide two different mass splittings. 

The loop corrections, shown in fig. |1|, from slepton-lepton and squark-quark loops to 
each element of the mass matrix are given by ||39[| , 



[m. 



1 / M'-^'' M*'^'* 

gm jloop = J2 Kkp^mpkmP slu 24>f^ lu — + 3 ^ Kkp^rapk'^T 20^^^ lu — ^ 



Here the angles (p^f^ and 0^^-' are given in terms of the slepton and squark mass eigenstates 
for a given flavour (p) Mi 2 and Mgi^q2 as 

,i„ 2< In = m? f-pH^-^Ln \^]]= mT' { T^fi^D 1 (8) 




where rafq are lepton and quark masses of flavour p; X = A + ^ tan P , which we take 
to be universal, and provides a mixing term between the left- and right-handed squarks and 
sleptons. The function / is defined to be 



"^This can be done order by order in the loop expansion when one appropriately defines the mass 
matrices of the Higgs sector . 
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We consider here that we are in the down-quark mass eigenstate basis, and that the A, A' 
couphngs have been redefined in terms of the couphngs appearing in the superpotential and 
of the corresponding Cabibbo-Kobayashi-Maskawa matrix elements; however, for simphcity 
we do not introduce additional notation. 

In this model, in order to simplify, we drop the slepton-fiavour dependence from 
i fi^e'^) ^^'i j^,j(p) i fi^'f^) and consider them to be universal in the slepton and squark 
sector respectively. We are then led to 

(m,m)ioop = 77^^ ( H KkpXmpkmf^m'f'' + 3^^^ Kkp^mpk^f^ ] • (11) 
\ 2 k,p q2 k,p J 

We use the following parametrization for the one-loop mass matrix written in the same basis 
as the tree-level mass matrix: 



loop 



/ mil mi3 \ 

mi2 17122 17123 ■ (12) 

\ mis 1^23 / 

The loop mass matrix elements niij are given in appendix 0, where we have kept only the 
relevant contributions after employing the mass hierarchy selection from the charged-lepton 
and down-quark sector. 

Finally, the total mass matrix we consider is 

M, = M'J'' + M^°°^. (13) 

An important remark is that there are additional loop corrections that can generate 
neutrino masses through the sneutrino mass splitting in these models f\ An additional 
gaugino-sneutrino (Higgs) loop appears where the R-parity violation is now contained in 
the mass-mixing term between the sneutrino and antisneutrino. Naively, it can be shown 
that this effect is small when the soft mass vector B = B^^ and the vacuum expectation 
value vector of the down-type scalar fields Va are aligned [|62| , |63[ . However, a full treatment of 



this effect requires a one-loop calculation of the scalar mass matrices (Higgses and sleptons), 
which is beyond the scope of this paper. 

A. Toy model 

In order to have an idea of how to constrain the whole space of parameters given by the 
different A and A' involved in the loop contribution to the mass and the /ie,^,r appearing in 



^We especially thank S. Davidson for discussions on this point. 
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sleptons 



leptons 



squarks 



quarks 



FIG. 1. One-loop Feynman diagrams contributing to the neutrino masses. 

the tree-level mass with the neutrino data, we put ourselves in the naive picture where all 
the A couplings are equal and also the A' couplings. In fact we are naively assuming in eq. 
(|ll|) (see appendix 0) that A133 = A233 = A333 and A'133 = A233 = A'333. 

We thus have two universal trilinear R-parity couplings, A and A', and three bilinear cou- 
plings, fie, /i^ and fir- In this toy model, our mass matrix can be written as follows 



/ Ki + K2 + al Ki + K2 + a^a^j, K2 + a^a^ \ 
Ki + K2 + OeU^ Ki + K2 + af^ K2 + a^Or 

\ K2 + Oettr K2 + a^ttr K2 + J 



(14) 



where 



Generically this matrix has three non-zero eigenmasses. Under several different approxi- 
mations (i.e. Ki = 0, or K2 = 0, or = Og), this mass matrix yields only two non-zero 
eigenvalues 

[/* X A< X f/ = diag{0,m2,m3}. (16) 

The different physical motivations (or limits) for these approximations correspond to the 
cases where: 

• A = or /(x^) = f{xq) and M2 = gives Ki = due to the hierarchy nib ^ "^r- 

• A' = or, if we take the limit in which the squarks decouple, K2 = 0. 
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• assume some hierarchy between the bihnear terms in the Lagrangian, such that ~ 



V. RESULTS 



A. General results 
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FIG. 2. Am^ in units of eV^ as a function of sin 2Q. The atmospheric and different solar 
solutions that are presented result from a combined fit with CHOOZ and SuperK constraints. 
The atmospheric solution is represented by the parameters (Amf3, sin^2023)- The MSW-LMA 
and MSW-SMA solar solutions are represented by the parameters (Amf2 or Am23 and sin^20i2) 
while the vacuum solution is represented by the parameters (Amf2, sin^ 2^12). 

We have performed a general scan of parameter space made up by the six matrix elements 
that appear in eq. (H). In this section we will simply present the results, and our discussion 
of them is given in the next section. We present in fig. ^ the allowed points on the (Am^, 
sin^ Iff) plane, which satisfy the combined constraints (CHOOZ, SuperK and solar data) for 
the possible solutions of the solar and atmospheric anomalies. Allowing only combined so- 
lutions of solar, SuperK and CHOOZ data restricts the available region of parameter space. 
We also see in fig. |^ that all three possible solutions for the solar anomaly, CHOOZ and 
atmospheric SuperK bounds can be simultaneously satisfied. After several runs, we choose 
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the interval of variation of the matrix elements in such a way that, combined to CHOOZ, 
SuperK constraints and one (or all) of the solar solutions is (are) possible. For example, 
in order to have a combined fit of the vacuum solution for solar anomaly, together with 
CHOOZ and SuperK constraints then the matrix elements have been varied in the interval 
-0.5 eV < Mij < 0.5 eV. 

We also note that the density of points (solutions) for Vacuum + CHOOZ + SuperK is 
smaller than the one for MSW-SMA + CHOOZ + SuperK which in turn is much smaller 
than the one corresponding to MSW-LMA + CHOOZ + SuperK. 

Tables HI, IV and V, which are given in the next subsections, present the allowed ranges 
of the input parameters (matrix elements) that satisfy the different solar oscillation solutions 
together with the CHOOZ and the atmospheric SuperK constraints. 

To study the spectrum for each of the possible combined solutions we give samples of 
the different mass spectra for different values of the matrix elements Mij. For illustration, 
we have present our plots for the spectrum as a function of the matrix element M13, which 
varies when we scan over parameter space. We also plot the effective mass relevant for 
neutrinoless double beta decay and the sum of the mass eigenvalues. Our results for the 
latter two quantities show that for the specific range that we have chosen to satisfy the 
combined fit of solar, CHOOZ and atmospheric constraints, the present bounds are always 
satisfied. In fact, for Major ana neutrinos, a shghtly stronger bound on the effective mass 
can be placed by our analysis and will be given below in each subsection. This bound is 
valid when we vary the mass matrix elements in the ranges that are specified below for each 
subcase. 
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FIG. 3. Am^ in units of eV^ function of sin 29. The atmospheric and solar solutions 
that are presented result from a combined fit of MSW-LMA, CHOOZ and SuperK constraints. 
The atmospheric solutions are represented with the parameters (Am^g, sin^2023) while the 
MSW-LMA solar solutions are presented with the parameters (Amf2 or Am23 and sin^20i2). 

1. MSW-LMA 

For a given model, imposing MSW-LMA, atmospheric SuperK and CHOOZ constraints 
simultaneously determines the allowed ranges for the input parameters Mn (eV), M22 (eV), 
M12 (eV), M23 (eV), M33 (eV), M13 in (eV), summarized in table IIL In order to obtain these 
results we have initially allowed all six matrix elements to vary in the interval [—1.0, +1.0] 
eV. 

Figure ^ shows the available regions in parameter space from the combined constraint of 
SuperK data together with CHOOZ and the MSW-LMA solution. 

It is clear from fig. |^ that, with this scan of parameter space, we can cover most of the 
allowed region that satisfies the combined constraints mentioned above. 

In figs. ^ and |b we plot the results of the individual mass eigenvalues and the effective 
mass (only valid for Majorana masses) together with the sum of the eigenmasses, which is 
the relevant quantity for HDM when we impose the constraints simultaneously. We can see 
from fig. ^ that the hierarchical spectrum is preferred for small values of our parameter 
M13. 

Using the results presented in fig. Qd, for the entire range of variation of the matrix element 
Mi3, a strong bound can be placed on the neutrinoless double beta decay effective mass, i.e. 
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(a) lb) 
FIG. 4. The spectrum presented results from a combined fit with CHOOZ and SuperK and 

MSW-LMA solar constraints. Mass spectrum (a) which presents the three possible mass eigenvalues 

(eV). The effective mass and the sum of eigenmasses are given in (b). 

for Mi3 G [—4, +4] x 10~^ eV (while the other matrix elements vary in the intervals given 
in table III), ;$ 5 x 10"^ eV. 



2. MSW-SMA 



Imposing MSW-SMA, SuperK and CHOOZ constraints simultaneously gives the ranges 
for the input parameters summarized in table IV. As in the previous section, this illustrative 
table is valid for a given model in which the mass matrix elements vary initially, before 
applying the solar, atmospheric and CHOOZ constraints, in the range [—1.0, +1.0] eV. 

Figure §a gives the combined constraint coming from SuperK data together with the 
CHOOZ bound and the MSW-SMA solution. We have noticed that with our scan of the 
parameter space the resulting density of points in fig. ^ is much smaller than in the case 
of MSW-LMA (fig. I). 

We show in fig. ^b the results of the individual mass eigenvalues, the effective mass 
and the values of Smj when we impose the constraints. Here once again there is a bound 
rUes ~ few X 10~^ eV for M13 e [—4.0, +4.0] x 10~^ eV (the other matrix elements take 
values in the ranges given in table IV). 



14 



Inputs 


Minimum 


Maximum 


Mii(10-2eV) 


-2.7 


7.3 


Mi2(10-2eV) 


-4.3 


4.3 


Mi3(10-2eV) 


-3.9 


3.9 


M22(10-2eV) 


-1.1 


8.3 


M23(10-2eV) 


-7.4 


7.4 


M33(10-2eV) 


-5.0 


8.9 



TABLE III. Allowed ranges of the input parameters that satisfy MSW-LMA, SuperK and 
CHOOZ constraints simultaneously. They also fulfil neutrinoless double beta decay and cosmolog- 
ical constraints. This illustration is for a model where the allowed inputs vary initially in a given 
interval, i.e., Mij G [-1.0, +1.0] eV . 



3. Vacuum Solar Solution 

In this case, there are few solutions that satisfy simultaneously all of the constraints. One 
can see in fig. |^ how the allowed points on the (Am^, sin^2^) plane are reduced compared 
to the density of points in the cases of MSW-LMA (fig. |) and MSW-SMA (fig. |a). 
The allowed ranges for the input parameters, after imposing the vacuum oscillation solution, 
SuperK and CHOOZ constraints simultaneously, are summarized in table V. This illustrative 
table is valid for a given model in which the mass matrix elements vary initially, before 
applying the solar, atmospheric and CHOOZ constraints, in the range [—0.5, +0.5] eV. 
In fig. where the spectrum of eigenmasses is represented together with the effective 
mass meff and the sum over the eigenmasses, a salient feature is the presence of two nearly 
degenerate eigenvalues (mi ~ m2). Here we also have a strong bound on the effective mass 
mefT ;$ fewlO"^ eV as we vary M13 in [—4.0,4.0] x 10~^ eV (the other matrix elements take 
values in the ranges given in table V). 
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Inputs 


Minimum 


Maximum 


Mii(10-2eV) 


-4.0 


2.9 


Mi2(10-2eV) 


-4.0 


4.0 


ivi 13 V ^ c V I 


—4 n 


4 n 


M22(10-2eV) 


-1.3 


8.3 


M23(10-2eV) 


-8.0 


8.0 


M33(10-2eV) 


-4.0 


5.9 



TABLE IV. Allowed ranges of the input parameters that satisfy MSW-SMA, SuperK and 
CHOOZ constraints simultaneously. They also fulfil neutrinoless double beta decay and cosmo- 
logical constraints. This illustration is for a model where the allowed inputs vary initially in the 
interval [-1.0, +1.0] eV . 
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(a) (b) 
FIG. 5. In (a) we represent Am^ in units of eV^ as a function of sin^ 29. The atmospheric 

solutions are represented with the parameters {Aml^, sin^ 2^23) and the MSW-SMA solar solutions 

with the parameters (Am^g oi' ^"^23 sin^2^i2) . In (b), the mass spectrum is represented 

together with rries and the sum of the eigenvalues (eV) . The atmospheric and solar solutions that 

are presented result from a combined fit with CHOOZ and SuperK constraints. 
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FIG. 6. In (a) we represent Am^ in units of eV^ as function of sin^ 29. The square and the 

triangle refer to the atmospheric solution {Amf^, sin^ 2^23) and the vacuum solar solution (Amfg, 

sin^20i2), respectively. In (b), the mass spectrum is represented together with meg and the sum 



of the eigenvalues (eV). 



Inputs 


Minimum 


Maximum 


Mii(10-2eV) 


-0.7 


1.3 


Mi2(10-2eV) 


-4.0 


4.0 


Mi3(10-2eV) 


-4.0 


4.0 


M22(10"2gY) 


+0.98 


4.1 


M23(10-2eV) 


-3.5 


4.1 


M33(10-2eV) 


+0.98 


4.1 



TABLE V. Allowed ranges of the input parameters that satisfy solar Vacuum, SuperK and 
CHOOZ constraints simultaneously. They also fulfil neutrinoless double beta decay and cosmo- 
logical constraints. This illustration is for a model where the allowed inputs vary initially in the 
interval [-0.5, +0.5] eV . 
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B. Bounds on R-parity-violating couplings from neutrino data 

Our approach to the numerical analysis is to use low energy input without any reference 
to high-scale physics. 

The loop contributions are proportional to different R-parity-violating parameters from 
those which control the tree-level terms. Three different possibilities can be considered: 

- the tree-level contributions are much larger than the loop contributions, 

- the tree-level contributions are of the same order as the loop contributions, 

- the tree-level contributions are much smaller than the loop contributions. 

Applying the results of our general scan of parameter space, we are in fact including all 
of the above-mentioned cases in our analysis. As each combination of constraints defines 
a different allowed range of parameters we will present below the corresponding bounds on 
R-parity-violating couplings separately for each case. 

1. MSW-LMA 



Imposing MSW-LMA, SuperK and CHOOZ constraints simultaneously defines the ranges 
for the input parameters Ki (eV) (which is proportional to A^), K2 (eV) (which is propor- 
tional to A'^), fte (eV^/^), (eV^/^), (eV^/^), which we have summarized in table VI. 



Inputs 


Minimum 


Maximum 


i^i(eV) 


0.0 


3.0 X 10-2 


K2ieV) 


0.0 


2.7 X 10-2 


ae (eVV2) 


-0.114 


0.154 


(eVV2) 


-0.280 


0.28 


ar (eVV2) 


-0.238 


0.238 



TABLE VI. Allowed ranges of input parameters that satisfy MSW-LMA, SuperK and CHOOZ 
constraints simultaneously. They also fulfil neutrinoless beta decay and cosmological constraints. 



Figure |^ shows the regions that satisfy the combined constraint coming from SuperK 
data together with CHOOZ and the MSW-LMA solution. 

We plot in fig. ||a the results of the individual mass eigenvalues, and in fig. the 
effective mass meg together with the sum over the eigenmasses when we impose all the 
constraints simultaneously. For illustration, we show how the mass spectrum and the ef- 
fective mass vary as a function of x = |ae(eV^^2)|, One can read from fig. ^ that in the 
range of parameters x G [0,0.1], we can put a strong bound on the effective mass, that is 
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FIG. 7. Am^ in units of eV'^ as a function of sin^ 26. The square and the circle refer to the 
atmospheric solution {Aml^, sin^2023) and to the MSW-LMA solar solution (Amf2, sin^20i2), 
respectively. 

meff ^ 5 X 10^^ eV. We can also determine, in the same region of the parameter space, a 
strong limit on mi + m2 + < 0.2 eV. 



Applying our results from table VI to the specific mass matrix of eq. (|T^) we find that 
the following bounds on the Rp-violating parameters exist for each solar solution combined 
with SuperKamiokande and CHOOZ data constraints. For the trilinear couplings we have 



A' < X (3.0 X 10-2 



(17) 



y < . . X [2.7 X 10' 



3Xf{x,)ml'^'^- ^^^^ 
Using the results of the appendix M we can also give a bound on the bilinear couphngs. 



1^1 < 



4detM 



m| cos2 /3(Mi + tan2 ewM2) 

4detM 
m| cos2 p{Mi + tan^ 9wM2) 

AdetM 
m| cos2 (3{Mi + tan^ ^v^^Ma) 



X 



0.154 eV-^/^) 
X (0.28 eV-^/2 
(0.238 eV-^/2 



(19) 
(20) 
(21) 
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FIG. 8. The mass spectrum (eV) is represented in (a). The effective mass and the sum of 

eigenmasses (eV) are represented in (b). This is when the MSW-LMA solution, the atmospheric 

one and CHOOZ constraint are imposed simultaneously. 

2. MSW-SMA 



Imposing MSW-SMA, SuperK and CHOOZ constraints simultaneously determines the 
ranges for the input parameters given in table VII. 

Figure ^ gives the combined constraint coming from SuperK data together with CHOOZ 
and the MSW-SMA solution. As mentioned previously for our analysis of a generic mass 
matrix, the density of points is much smaller than in the case of the MSW-LMA solution 

(fig- 0)- 

As in the previous subsection by applying our results from table VII to the specific mass 
matrix of eq. ([I^ ) we find the following bounds on the Rp-violating parameters for each 
solar solution combined with SuperKamiokande and CHOOZ data constraints: 



Xf{xi)ml 
3Xf{x,)ml 



X 3.0 X 10 



3 



4detM 



m| cos^ /?(Mi + tan2 ewM2) 



(3.0 X 10 
(4.2 X 10 



eV-1/2 



(22) 

(23) 
(24) 
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Inputs 


Minimum 


Maximum 


i^i(eV) 


0.0 


3.0 X 10"'^ 




0.0 


3.0 X 10"^ 


ae (eVV2) 


-4.2 X 10-2 


4.2 X 10-2 




-0.28 


0.28 




-0.21 


0.21 



TABLE VII. Allowed ranges of the input parameters that satisfy MSW-SMA, SupcrK and 
CHOOZ constraints simultaneously. They also fulfil neutrinoless beta decay and cosmological 
constraints. 



< ^15777 X (0.28 , (25) 

n 4detM / ,/,,\2 , ^ 
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E le-04 



Atmospheric 



sin 2e 

FIG. 9. Am^ in units of eV^ as a function of sin^ 29. The square and the circle refer to 
the atmospheric sohition {Aml^, sin^2023) and the MSW-SMA solar solution (Am^2) sin^20i2), 
respectively. 

3. Vacuum 

In our present application, we find that imposing vacuum oscillations, atmospheric and 
CHOOZ constraints simultaneously or vacuum oscillation and CHOOZ constraints simul- 
taneously offers no solution of the solar and atmospheric anomalies. We have determined 
which are the allowed ranges for the input parameters by imposing only vacuum solution 
constraints; these are summarized in table VIII. 



Inputs 


Minimum 


Maximum 


i^i(eV) 


0.0 


6.0 X IQ-^ 


K2ieV) 


0.0 


1.20 X 10-2 


Oe (eVV2) 


-0.35 


0.35 


(eVV2) 


-0.35 


0.35 


ar (eVV2) 


-0.35 


0.35 



TABLE VIII. Allowed ranges for input parameters that satisfy Vacuum oscillation constraints. 
They also fulfil neutrinoless beta decay and cosmological constraints. 



We determine the bounds on the Rp-violating parameters from the results from table 
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VIII. These are 



Xf{xi)ml ^ ^ 



A < , p X (1.20 X 10-^1 

ml cos2 /?(Mi + tan2 OwM^) ^ ^ ^ ' 
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4- Limits on A, A' and jii 



Given our results for the bounds on the Rp-violating couphngs, the analysis in terms 
of Rp-conserving SUSY parameters is straightforward. Our general results can be analysed 
by considering the SUSY particle spectrum in detail; we will present these results in a 



forthcoming paper ||6J]. A common assumption used to place bounds in this model is to 
take all Rp-conserving mass parameters to be of the same order, M^usy For the particular 
case where Msusy = 100 GeV, /(x) 1 and tan/3 = 2, the bounds are 

• MSW-LMA + CHOOZ + SuperK 

< 1.5 X 10^^ , A'^ < 7.02 X 10"^ , 
f^e < 1.9 X 10"^ GeV, fXf, < 3.4 x 10"^ GeV, fir < 2.9 x 10'^ GeV . (30) 

• MSW-SMA + CHOOZ + SuperK 

A^ < 1.5 X 10-*^ , A'2 < 7.8 X 10-1° , 

He < 5.2 X 10"^ GeV, /i^, < 3.4 x 10"^ GeV, fir < 2.6 x 10"^ GeV . (31) 

• Vacuum 

A^ < 3.0 X 10"^ , A'2 < 3.1 X lO"'^ , 

fii < 4.3 X 10"^ GeV, i = e,fi,T . (32) 



C. Dependence of the (Am^, sin^ 26*) solar and atmospheric parameters on the 

Rp-violating parameters 

We now plot in figs. 0, |TT] and the different solar and atmospheric parameters as 
functions of the Rp-violating parameters Ki oc A^, K2 oc A'^, ae,^,r and also = al+af^+a^. 

In figs. and |ll], where the constraints from CHOOZ, SuperK and the MSW-LMA 
solution are satisfied simultaneously, we have plotted the solar parameters Am^2 sin^ 26*12 
(fig. |10D and the atmospheric variables Am^g and sin^2^23 (fig- |ll])- Figure |12| illustrates 



the dependence of the solar and atmospheric parameters on the sum of the squares of the 
Rp-violating bilinear terms for the combined fits of SuperK, CHOOZ and MSW-LMA (figs. 
12|a and |l2|c) as well as for Super, CHOOZ and MSW-SMA (figs. |T2|b and |l2|d). Here 



one can observe how the region of parameter space from the fit that satisfies MSW-LMA, 
CHOOZ and SuperK simultaneously is considerably larger than the available region where 
we fit the MSW-SMA, CHOOZ and SuperK atmospheric constraints. 
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VI. DISCUSSIONS AND CONCLUSIONS 



Our results from the general analysis show that there are many regions in parameter 
space that can accommodate simultaneously the MSW-LMA and atmospheric oscillation 
solutions together with the CHOOZ constraint, see figs. 2 and 3. In comparison, we see 
from fig. ^ that not all of the allowed region in the (Am^, sin^ 2^) plane is covered when 
MSW-SMA and the atmospheric oscillations solutions together with the CHOOZ constraints 
are imposed. Moreover, the region of parameter space for which the atmospheric, CHOOZ 
and vacuum oscillation constraints are satisfied has very few solutions. In contrast, we 
find that for the specific case of the MSSM with R-parity violation there are no regions in 
parameter which can accommodate simultaneously the vacuum solar solution and all other 
constraints. 

As far as the possible spectra are concerned, we clearly see that for the combined MSW- 
SMA, CHOOZ and SuperK constraints both the hierarchical and pseudo-Dirac spectra are 
more commonly present. There are also a few values of parameter space that present a de- 
generate spectrum. For a combination of MSW-LMA, CHOOZ and SuperK constraints, the 
three possible spectra are very common. As expected, for the combined vacuum oscillation 
solution, CHOOZ and SuperK constraints only pseudo-Dirac spectra are present. In refer- 
ence to figs. 4, 5b and 6b our conclusions hold and the bounds mentioned in the previous 
section are still maintained when we study the spectrum in terms of the other parameters 
Mij ^ Mil. 

It is clear that increasing the range of variation of our parameters changes the values of 
the eigenvalues of the mass matrix. In fact, larger values of the inputs make it impossible 
to satisfy the vacuum oscillation constraints; eventually the same will occur with the MSW- 
SMA and MSW-LMA solutions. It is really the combination of constraints that solve both 
atmospheric and solar anomalies which limit the allowed ranges. An important conclusion 
is that, after obtaining a simultaneous solution for the solar and atmospheric anomalies and 
the CHOOZ constraint, the effective mass constraint and the cosmological constraint are 
always satisfied. We have also indicated previously when a stronger bound can be placed 
on these quantities. 

Some of the previous work regarding neutrino masses in this specific model we consider 
have tried to fit only one of the Am^ to either the atmospheric or one of the solar solutions. 
From our analysis we have seen that fitting only one of the neutrino anomalies is straight- 
forward. What is much harder to obtain is a region in parameter space that satisfy the 
combined fits. Other authors have also studied the possibility of fitting both atmospheric 
and solar anomalies, although their approaches are somewhat different from ours. They con- 
sider either supergravity-inspired models or exclusively bilinear contributions, or exclusively 



25 



trilinear contributions. 

The usual bounds |]65| , |66| cited for the trihnear couphngs A133 and are placed when 



only one of these two R-parity- violating couplings are considered to be present in the theory. 
The bound is implemented by using the upper bound on the electron neutrino mass from 
the /^-decay spectrum. Having only one of the trilinears corresponds of course to having 
a single non-zero matrix element in eq. ([T^, given by Ki or K2 depending on whether the 
trilinear coupling is A133 or X[^^, respectively. We see that in the special limits we have 
considered on R-parity-conserving SUSY parameters our bounds from eq. ( ^Of ) agree with 
the previous bounds. The overall conclusion is that the bounds are very strong on the 
R-parity violating parameters. We have shown that the simultaneous inclusion of many of 
these couplings can constrain them much more than other bounds arising from other low- 
energy or high-energy processes p3|,|5B[]. The overall conclusion is that the bounds are very 



strong on the R-parity violating parameters. 

Figures 10-12 illustrate the different dependences of the solar and atmospheric parameters 
on the R-parity couplings. It is clear that for large values of the parameters there are no 
solutions. The discreteness of the allowed values of the parameters is due in part to the step 
we use in our scan of the parameter space. We also see, of course, that the dependence on 
a single parameter is complicated, since many solutions can be found for the same value of 
the variable. This is due to the large number of different parameters involved. 

The general global interpretation of the experimental data is that the SuperK seriously 
challenges MSW-SMA and Vacuum Oscillations solutions [ET^. This is in agreement with 



the analysis we have done. Indeed, the majority of the solar anomaly solutions we have 
that also satisfy both SuperK and CHOOZ are MSW-LMA. The future experiments SNO, 
BOREXINO, ICARUS, HERON with solar neutrinos as well as KAMLAND, K2K, MINOS, 
BooNE..., and many others, will provide additional decisive information to pin down the 
neutrino mysteries. 

Acknowledgments We thank S. Davidson for very interesting and useful discussions. 
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FIG. 10. The solar parameters Am|2 ^-'^d sin^ 2^12 as functions of the different R-parity pa- 
rameters (a), af_i/oir (b), ETi ( oc A^) (c) and K2 (oc A'^) (d), in the case of a combined fit with 
MSW-LMA, SuperK and CHOOZ constraints. 
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FIG. 11. The atmospheric parameters Amf^ and sin^ 2^23 as functions of the different Rp 

parameters Oe (a), an/a-j- (b), Ki (oc A^) (c) and K2 (oc A'^) (d), in the case of a combined fit with 



MSW-LMA, SuperK and CHOOZ constraints. 
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Charged and Neutral Fermion Mass Matrices 



APPENDIX A: CHARGED COLOUR-SINGLET MASS MATRIX 



Written in tlie basis \E' 



2 ,/,2 ,/,2 ,/,2 ^ ^Yie matrix 



M 



tree 








g2V sin f3 

g2V cos 13 

//I mi 




V 



/U2 
/"3 



\ 




ma 
ma / 



(Al) 



can be diagonalized explicitly to obtain the charginos and charged-lepton tree-level masses. 
However, our analysis of the neutrino spectrum requires the ^Uj's to be very small, thus 
their effect on the eigenvalues and corresponding mass eigenstates will be tiny. It is a good 
approximation for small R-parity violation to neglect the RPV contribution to the masses 
in the charged leptons and charginos 



APPENDIX B: TREE-LEVEL NEUTRAL COLOUR-SINGLET MASS MATRIX 



'■N 




V 



Ml 










giD sin /3 


g2V sin (3 


V2 


V2 


giv cos 13 


31?; sin (3 


V2 


V2 





















giv sin (3 
g2V sin (3 



-1^3 



31 V cos (3 
32 f COS (3 











-fll 










-fJ'2 








\ 






; 



(Bl) 



In the limit where the small, the effective neutrino mass matrix has the "see-saw" 

structure 



M 



tree 



(B2) 



The assumption we make so as to include the one-loop corrections is that only the corrections 
to m° need be considered, thus giving the one-loop-corrected expression for the 3x3 neutrino 



mass matrix 42 , 5C 



(B3) 



30 



where M^°°p = ml + mJ,°°P. The expression for m[^°P is given in the next subsection. The 
elements of the matrix M^'^^ = ^M~^^'^, which is parametrized in terms of the a^'s of eq. 
(^) are given by 

Mr = d^''''Vf\T''^mvl (B4) 

as we are working in the basis where the slepton vacuum expectation value is zero, vi = 
Vd = f cos/3. The basis-invariant expression corresponds to substituting fii^jvf with < 
Vi > —ij,iVi){fi{vj) — iJ-jVi), where (vi) is the vacuum expectation value of the slepton fields 
LR 



APPENDIX C: LOOP-MASS MATRIX 

Expliciting the indices in eq. (|^ leads to the following symmetric matrix elements of 
which we keep only the relevant contributions due to the mass hierarchy in the charged- 
lepton sector nie -C <^ rrir and in the down-quark sector nid -C <ti nib. 

^11^ = TFZ^ ( (^?33^r) + S^TT^^ (A'133'^b + A'i23Ai32"^."^fe) ) (CI) 



loop 



m 



loop 



{^Mf (^133^r) + (A'lSS'^fe + A'i23Al32"^."^fe) j 

167r2 {^Mf ('^133A233rn^) + 3"^^ (A'i33A233"^6 + (^'132^223 + K23^232)msmb)^ (C2) 



167r2 
X 



loop _ 3 X /(x, 



167r2 Mf 



A'l33A333"^6 (C3) 



X 



167r2 {^^Mf ('^233'^r) + 3"^^ (A233"^b + A223A232m,mfe)') (C4) 



loop S X f{xg) f f 2 (r''^\ 

"^23 — j^q2 ^233^333"^?) \^^) 



loop 3 X f{Xq) ,2 2 fr^r\ 



For a general discussion of basis-independent parametrizations of R-parity violation, see refs. 



31 



REFERENCES 



[1] Y. Fukuda et al., Super-Kamiokande Collaboration, Phys. Lett. B 433 (1998) 9; Phys. 
Lett. B 436 (1998) 33; Phys.Rev.Lett. 81 (1998) 1562. 

[2] R. Becker-Szendy et al., 1MB Collaboration, Nucl. Phys. B (Proc. Suppl.) 38 (1995) 
331. 

[3] W.W.M. Allison et al. Soudan-2 Collaboration, Phys. Lett. B391 (1997) 491, Phys. 
Lett. B449 (1999) 137. 

[4] Y. Fukuda et al, Kamiokande Collaboration, Phys. Lett. B 335 (1994) 237. 

[5] R. Davis et al, Phys.Rev.Lett. 21 (1968) 1205; B.T. Cleveland et al., Astrophys. J. 496 
(1998) 505. 

[6] W. Hampel et al., Phys. Lett. B 388 (1996) 384. 

[7] D.N. Abdurashitov et al, Phys. Rev. Lett. B 77 (1996) 4708; |astro-ph/ 9907131 



K. S. Hirata et al., Kamiokande Collaboration, Phys. Rev. Lett. 77 (1996) 1683. 
[9] Y. Fukuda et al., Super-Kamiokande Collaboration, Phys.Rev.Lett. 81 (1998) 1158. 
[10] M. Ambrosio et al., MACRO Collaboration, Phys. Lett. B 434 (1998) 451. 
[11] C. Anthanassopoulos et al., LSND Collaboration, Phys. Rev. Lett. B 81 (1998) 1774. 
[12] S. Coleman and S.L. Glashow, Phys. Rev. D59 116008 and references therin. 
[13] A. Halprin and H.B. Kim. KIAS-P99028, preprint [hep-ph/9905301| and references therin. 
[14] S. Pakvasa, |hep-ph/990542B. 



[15] G.L. Fogli, E. Lisi, A. Marrone and G. Scioscia, Phys. Rev. D60 053006 and Phys. Rev. 
D59 117303. 

[16] B. Zeitnitz et al.,KARMEN Collaboration Prog.Part.Nucl.Phys. 40 (1998) 169. 

[17] see for example and refs. therein: R.Barbieri, J. Ellis and M.K. Gaillard, Phys. Lett. B90 
(1980) 249; M. Gell-Mann, P. Ramond and R. Slansky, Proceedings of the Supergravity 
Stony Brook Workshop, New York 1979, eds. P. Van Nieuwenhuizen and D. Freedman; 
T. Yanagida, Proceedinds of the Workshop on Unified Theories and Baryon Number in 
the Universe, Tsukuba, Japan 1979, ed.s A. Sawada and A. Sugamoto; S. Dimopoulos, 
L.J. Hall and S. Raby, Phys. Rev. Lett. 68 (1992) 1984; G.K. Leontaris et al, Phys. 
Rev. D53 (1996) 6381; S. Lola and J. D. Vergados, Prog. Part. Nucl. Phys.40 (1998) 
71; G. Aharelli and F. Feruglio, |tiep-ph/9905536| and |hep^h/9907532| ; J.A. Casas et al. 



32 



hep-ph/9905381| and |hep-ph/990628ll ; M. Carena et al, |hep-ph/9906362| ; P. Langacker, 
Nucl. Phys. Proc. Suppl. 77 (1999) 241; Z. Berezhiani and A. Rossi, JHEP 9903 (1999) 
002. 

[18] R. Barbieri, L.J. Hall, D. Smith, A. Strumia and N. Weiner, JHEP 9812 (1998) 017. 
[19] G. Gelmini and E. Roulet, Rept. Prog. Phys. 58 (1995) 1207. 

[20] B. Kayser, F. Gibrat-Bebu and F. Perrier, The Physics of Massive Neutrinos, Wolrd 
Scientific, Singapore 1989. 



[21] S.M. Bilenkii, C. Giunti, W. Grimus, |hep-ph/9812360 
[22] J.W.F. Valle, ^ep-ph/98092"^ . 

[23] M. ApoUonio et al., CHOOZ Collaboration, Phys. lett. B 420 (1998) 397. 
[24] H. Georgi and S.L. Glashow, |hep-ph/ 980829^ . 

[25] M. Giinther et al., Phys. Rev. D55 (1997) 54, Phys. lett. B 407 (1997) 219. 
[26] L. Baudis et al., Phys. Lett. B 407 (1997) 219. 

[27] R.A.C. Croft, W. Hu and R. Dave, |astro-ph / 9903335| , and see references therein. 

[28] V. M. Lobashev, to be published in the proceedings of 17th International Workshop on 
Weak Interactions and Neutrinos (WIN 99), Cape Town, South Africa, 24-30, Jan 1999; 
Prog.Part.Nucl.Phys.40 (1998) 337. 

[29] H. Barth, Prog.Part. Nucl. Phys. 40 (1998) 353, J. Bonn, to be published in the proceed- 
ings of 17th International Workshop on Weak Interactions and Neutrinos (WIN 99), 
Cape Town, South Africa, 24-30, Jan 1999. 

[30] G. Gelmini, to be published in the proceedings of 17th International Workshop on Weak 
Interactions and Neutrinos (WIN 99), Cape Town, South Africa, 24-30, Jan 1999, |hep- 
I ph/ 9904369 . 

[31] G.L. Fogli, E. Lisi, A. Marrone and G. Scioscia, Phys. Rev. D59 33001. 

[32] V. Barger, T.J. Weiler and K. Whisnat, Phys. Lett. B440 (1998)1. 

[33] M. C. Gonzalez-Garcia, H. Nunokawa, O.L. G. Peres and J.W.F. Valle, Nucl. Phys. 
B543 (1999) 3. 

[34] J.N. Bahcall, P.I. Krastev and A.Yu. Smirnov, Phys. Rev. D58 96016. 
[35] N. Hata and P. Langacker, Phys. Rev. D56 (1997) 6117. 



33 



[36] S.T. Petcov, to be published in the proceedings of 17th International Workshop on 
Weak Interactions and Neutrinos (WIN 99), Cape Town, South Africa, 24-30, Jan 1999, 
|hep-ph/99TI7^ . 



[37; 
[38 
[39 
[40 
[41 
[42 
[43 
[44; 
[45 
[46 
[47 
[48 
[49; 
[50 
[51 
[52 
[53 
[54 
[55 
[56 
[57 
[58 
[59 



L.L. Chau and W.Y. Keung, Phys.Rev.Lett. 53 (1984) 1802. 
R. Godbole, R Roy and X. Tata, Nucl. Phys. B401 (1993) 67. 
Y. Grossman, H.E. Haber Phys.Rev. D59 093008. 
E.Nardi, Phys. Rev. D55 (1997) 5772. 

T. Banks, Y. Grossman, E. Nardi and Y. Nir, Phys. Rev. D52 (1995) 5319. 

M. Nowakowski and A. Pilaftsis, Nucl. Phys. B461 (1996) 19. 

P.M. Borzumati, Y. Grossman, E. Nardi and Y. Nir, Phys.Lett. B384 (1996) 123. 

R. Hempfling Nucl.Phys. B478 (1996) 3. 

B. de Carlos and P. White, Phys. Rev. D54 (1996) 3427. 

A.Y. Smirnov and F. Vissani, Nucl. Phys. B460 (1996) 37. 

H.P. Nilles and N. Polonsky, Nucl. Phys. B499 (1997) 33. 

D. E. Kaplan and A.E. Nelson, |hep-ph/9901254| . 

E. J. Chun and J.S. Lee, |hep-ph/98ll2UT. 



V. Bednyakov, A. Faessler and S. Kovalenko, Phys. Lett. B442 (1998) 203. 
J.C. Rom&o, |hep-ph/ 9907466 . 

S.Y. Choi, E.J. Chun, S.K. Kang and J.S. Lee, |hep-ph/99034H3 . 

M. Dress, S. Pakvasa, X. Tata and T. ter Veldhuis, Phys. Rev. D57 (1997) 5335. 

S. Rakshit, G. Bhattacharyya and A. Raychadhuri, Phys. Rev. D59 091701. 



G. Bhattacharyya, H.V. Klapdor-Kleingrothaus, H. Pas, |hep-ph/ 9907432 



E.J. Chun, S.K. Kang, C.W. Kim and J.S. Lee, Nucl. Phys. B544 (1999) 89. 
K. Choi, E.J. Chun and K. Hwang, |hep-ph/9811363 . 



O.C.W. Kong, Mod. Phys. Lett. A14 (1999) 903. 

A. Joshipura and S. Vempati, |hep-ph/9903435| . A. Joshipura, V. Ravindran and S. 
Vempati, Phys. Lett B451 (1999) 98; A. Joshipura and M. Nowakowshi, Phys. Rev. 



34 



D51 (1995) 2421. 
[60] A. Datta, B. Mukhopadhyaya and S. Roy, |hep-ph/9905549 



[61] H. Haber and G. Kane, Nucl. Phys. B272 (1986) 1. 
[62] S. Davidson and S. King, Phys. Lett. B445 (1998) 191. 

[63] Y. Grossman and H. Haber, Phys. Rev. Lett. 78 (1997) 3438; |hep-ph/99063TO| . 
[64] A. Abada and M. Losada, in preparation. 
[65] H. Dreiner, |hep-ph/97071H3. 



[66] G. Bhattacharyya, Nucl. Phys. Proc. Suppl. 52A (1997) 83 and |hep-ph/9709395 



[67] Chang Kee Jung, Talk given at International Europhysics Conference on High-Energy 
Physics (EPS-HEP 99), Tampere, Finland, 15-21 Jul 1999. 

[68] M. Bisset, Otto C.W. Kong, C. Macenasu and Lynne Orr, |hep-ph/9811"498| (1998). 

[69] S. Davidson and J. Elhs, Phys. Lett. B39 (1997) 210; Phys. Rev. D56 (1997) 4182. 



35 



